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1. Introduction 

One of the major issues in QFT is the search for non-perturbative results. In particular little is known 
about the phase structure (fixed points, critical exponents, etc.) of most physically interesting theories. Nor 
have intriguing questions about the renormalization flows, like the proper extension of the c-theorem [1], 
been completely understood [2]. 

One of the methods capable of handling such problems is the exact renormalization group (ERG, here- 
after). Originally developed by Wilson in his seminal articles in the early seventies [3] (see Ref. [4] for a 
classical review), it has recently attracted much attention. Although firstly used for studies of critical phe- 
nomena in condensed matter problems [5] its scope includes many other fields like particle theory, as it was 
demonstrated by Polchinski in an elegant paper where he proved the perturbative renormalizability of X(j) 4 
theory in a quite simple way [6]. Similar manipulations have led to quite interesting results regarding the 
study of different aspects of the perturbation expansion around a Gaussian fixed point and its associated 
diagramatic expansion [7,8]. However, the power of the ERG relies on the possibility of obtaining quantita- 
tive knowledge about the renormalization flows, in particular one piece of information which is probably the 
most valuable, due to its universality: the critical exponents [9]. 

The ERG approach is based on writing down a functional differential equation that expresses how the 
action changes when we integrate out high energy modes. This is the so-called ERG equation, and it is 
the cornerstone of the whole technique: with it and together with the most general action consistent with 
the symmetries of the model, the complete set of /3-functions can be computed, and from these the location 
of the fixed points and their exponents. However, because of practical reasons it is impossible to handle 
all possible operators which could be included into the action. One must choose a more selective criterion, 
rather than simply to be consistent with the demanded symmetries, that is, one must choose a reasonable 
truncation of the general expansion. Usual approximations attempt to restrict the space of interactions to a 
reasonable number of operators, e.g. replacing the effective action for a non-derivative effective potential or 
expanding the action in powers of the momentum. For instance, to study the Wilson fixed point of A</> 4 in 
three dimensions, the authors of Ref. [9] consider only arbitrary polynomials in the fields without derivatives, 
but the same type of problem has also been addressed by changing the ERG equation and/or by considering 
other types of truncations [4,10-19]. 

ERG methods have been used in more complicated cases, like phase transitions at finite temperature 
[20], theories with gauge interactions [8,21] and theories with fcrmions as fundamental particles. In the last 
case, it usually relies on a certain kind of bosonization, usually consisting in coupling fcrmionic bilinears with 
scalar fields, and then studying self-interactions of these scalar fields 1 [22]. We feel, nonetheless, that this 
approach is unsatisfactory. The reason is that one should learn how to deal with pure spinor theories without 
simplifying the problem to a scalar one. Moreover, there are some quite interesting phenomenological models 
described entirely with spin 1/2 fields, like the celebrated Fermi theory of weak interactions [25], models for 
resonance physics [26] based on extensions of the Nambu-Jona-Lasinio action [27], or even models explaining 

1 The main exceptions are Ref. [23] which deals with perturbative properties of the Gaussian fixed point and the 
series of Ref. [24] which lead to the proof of the existence of the Gross-Neveu model in 2 + e dimensions. We are 
seeking, however, non-perturbative quantitative information out of the ERG, rather than rigorous formal statements. 
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the symmetry-breaking sector of the electroweak theory, especially in connection with technicolor theories 2 
[29] . It is also compelling the curiosity that all fields in the Standard Model can be expressible as fermionic 
fields [30]. Surely, further extensions have to be dealt with afterwards, like scalar and spin 1/2 fields coupled 
together and spinor particles interacting through gauge fields. 

On the other hand, we remind the reader that fermions are not always easily manageable by non- 
perturbative methods, e.g. Lattice Field Theory. On the contrary, we will show that, once truncated, the 
ERG equation treats fermions and bosons very similarly, thus making possible a nearly immediate translation 
of knowledge from one case to the other. 

Our purpose is to work with a sample theory based solely on spinor fields and to develop a method of 
obtaining numerical non-perturbative information from it (e.g. some critical exponents). With this motiva- 
tion an ERG equation is derived, similar to the bosonic one by Polchinski [6], and applied to a particular 
model. We try to emphasize throughout the article that many of the peculiarities encountered are nothing 
more than the translation of their counterparts already found in previous papers on the ERG for bosonic 
theories. 

To begin with, one may choose an appropriate model, relatively simple and non-trivial. As usual, the 
two-dimensional world is a perfect site where to look for. Indeed, the two-dimensional Clifford algebra is the 
simplest one, generated by the well-known Pauli matrices. 

Moreover, the study of self-interacting fermionic theories in d — 2 can be traced back to the work of 
Thirring [31], where he proposed a massless model of a single Fermi field containing a quartic self-interaction. 
It can be solved exactly [32] and it presents some interesting features, most of them probably unexpected. It 
is perturbatively renormalizable, as power counting arguments suggest, but what is not expected from naive 
arguments is that it does not describe an interacting theory, but a trivial one [33]. 

The model becomes non-trivial when N species of fermions transforming under a global representation 
of the unitary group U(N) are considered. This is the Gross-Neveu model [34], which is asymptotically free 
and renormalizable, within perturbation theory and also within the 1/N expansion. However, none of these 
approximations is capable to find any non-trivial fixed point for d = 2. (Actually the 1/N expansion shows 
a non-trivial fixed point of order d — 2 in d dimensions). 

An interesting modification leads to the so-called chiral Gross-Neveu model [34], which is chosen to 
have the additional symmetry of the Ur(N) x Ul{N) chiral group. As in the previous case no fixed points, 
besides the Gaussian one, can be found within the 1/N approximation. Nevertheless, other non-perturbative 
techniques are available. Thus the quartic interaction of the chiral Gross-Neveu model can be expressed, 
after a Fierz transformation, as a current-current interaction and the latter allows an operator analysis of 
the model within the current algebra approach. Such a study is carried out in Ref. [35] where, exploiting 
conformal techniques, two critical curves in the space of couplings are found for which the theory is scale 
invariant. One of the lines corresponds to the abelian Thirring model, whereas the other one is truly non- 
trivial and does not pass through the origin. A very remarkable fact is that this result is exact and is 
not given by a zero of a /3-function, neither in perturbation theory nor in the large N expansion. For this 
continuum set of critical theories the value of the coupling constant, associated to the abelian degrees of 
freedom, is arbitrary while the coupling associated to the SU(N) degrees of freedom is fixed to be equal to 

2 Note, however, Ref. [28] 
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zero or 4ir/(N + 1). It is important to notice, however, that the fixed point is not unique and, as it happens 
in the Thirring model too, depends on an arbitrary parameter related to rcgularization ambiguities. 

More recently, using bosonization, current algebra and conformal techniques, other non-trivial fixed 
points in two-dimensional fermionic models were found [36] . However a lot of work has to be done in order 
to gain a clearer understanding of the complete phase diagram. 

Before ending this section let us describe shortly our results. We studied the chiral Gross-Neveu 
model with a Polchinski-type ERG equation and we projected the space of local interactions onto the 
106-dimensional subspace generated by terms with at most six fcrmions and three derivatives. 

Within the large N approximation we found two different non-trivial N — > oo limits, leading to qual- 
itatively different results. One of them leads to a continuous family of fixed points along the direction of 
the U(l) excitations, similar to that of Dashcn and Frishman. However, unlike their case, the anomalous 
dimension r\ vanishes at leading order in l/N. The critical exponents can be computed analytically and 
most of them coincides with the canonical values. However for a wide class of regularization schemes the 
most relevant critical exponent is non-trivial and takes the value Ai = 1.1231.... The other solution gives a 
non-trivial (but scheme dependent) anomalous dimension (rj = 1.0 ~ 1.4 in the range studied), but the fixed 
point is isolated. The relevant critical exponent in the same range is in this case Ai = 2.18 <~ 2.26. 

For finite N the fixed point equations can be solved only numerically. We found a plethora of solutions, 
most of them spurious due to the nature of the approximation. As we will discuss later, it is very difficult to 
discriminate between good and fictitious solutions. We took specially care of those solutions that for large 
N match the l/N fixed points. For both cases the solutions are strongly dependent on the scheme though 
the dependence can be mitigated using a minimum sensitivity criterion. In one case the solution is defined 
for any value of TV, with TV r\ growing asymptotically to 4.9, and Ai decreases to its l/N limit 1.123.... In the 
other case the solution disappears unexpectedly at N w 142 where it merges, as we will carefully describe, 
with another kind of solution. At the bifurcation point we have r\ = 1.88 and Ai = 5.8. 

Finally we analyzed the case N = 1. This case has to be treated separately because in absence of 
flavour the Fierz transformations impose additional constraints that reduce considerably the total number of 
independent operators. The results in this case are disappointing: the fixed point solution is isolated and not 
a continuous family as in the Thirring model. This property is fulfilled in the previous order approximation 
(terms with less than three derivatives) and it is unclear for us why it is lost at this order. For this case 
the values of the anomalous dimension has an extremely wide variation with the renormalization scheme 
(r,= 1-10). 

Let us remind the reader that in general the fixed-point actions are scheme dependent and they contain 
any possible operator, thus making it very difficult to interpret them directly. Nevertheless, any fixed point 
is characterized by that piece of information that is universal and, thus, physically relevant. This contains 
the number of relevant directions and their associate exponents. For these reasons, we have refrained from 
presenting in detail the fixed-point actions but, rather, we have concentrated in discussing the universal 
properties of our results. 

The article is organized as follows. In section 2 we derive an ERG equation for pure fermionic theories 
in any dimension. Sect. 3 is devoted to the construction of the action. The chiral Gross-Neveu model is 
defined through its symmetries and the truncation we use is explained. The calculation of the /3-functions, 
fixed points and the corresponding critical exponents is divided into sections 4 and 5, while the sixth one 
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contains a summary of results and some conclusions. Because of the unavoidable increasing complexity 
of the notation (we will be dealing with 106 operators) we have included an Appendix (A) where all our 
conventions are summarized. Although any new symbol is clearly defined when it appears, we thought it 
would be helpful to have all them condensed in a single page. We have also included one appendix (B) to 
write down the whole action and another one (C) to present the complete set of /3-functions. 



2. ERG equation for fermionic theories 

In this section we derive an ERG equation for a field theory written in terms of spinor quantities, on 
a Euclidean space of dimension d. Its role is to dictate the behaviour of the action as we integrate out 
modes, that is, how the action of our effective theory has to be modified when we vary the characteristic 
scale (cut-off) A, while keeping the S'-matrix elements invariant. More concretely, if we parametrize the 
renormalization flow with t = — In then the ERG equation will provide us with a sufficient condition to 
ensure Z = 0, A being a fixed scale and Z the generating functional for the connected Green's functions. 

We will consider, thus, a general theory whose action is expressible as a function of spinor fields only 3 
and artificially split it into 

S = Skin + Si n t ■ (2-1) 

Sint is an arbitrary polynomial in the fields and momenta (we will work always in momentum space) and 
Skin is a regulated version of the usual kinetic term, 



Skin= / iP_ p P^(p)iP P , (2.2) 
Jp 

where f stands for J d d p and Pa is the matrix 

Pa (p) - (2^) d ^A;$ , (2.3) 

with Ka (z) an analytic function over the whole finite complex plane that vanishes faster than any power-law 
when z — > +oo and is normalized to be Ka (0) = 1 [38] . 

In the following, we will consider expectation values defined by 

(X) = J I e - s f ^^/J-^a'W^+JJ-^a 1 Wx P +/« ; (2 4) 

where X is any operator, \ p , \ p are Grassmann sources, Qa (p) is another regulating function with analogous 
properties to Pa (jp) and, finally, /a is a c-number independent of the fields. 

With the above conventions, the starting point of the derivation is the observation that a functional 
integral of a total functional derivative vanishes, which leads us to 



tr (Pa^Pa) 6(0)) ((VVa" 1 - X-pQl 1 ) ^ (^V P - Ql'Xp)) , 



(2.5) 



3 We assume that Weinberg's conjecture [37] is valid: an arbitrarily general action leads to arbitrarily general 
S'-matrix elements and vice versa. 
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where the trace is over spinor indices. This is the counterpart of Eq. (1.9) of Ref. [38] and, as there, it can 
be used to identify the rate of change of the kinetic term, 



+ J ^-pP^PaQ^Xp + X-pQa'PaP^^p) + J (x-pQa'PaQ^Xp + tr (p^Pa) S (0)) . 



(2.6) 



On the other hand, by imposing that the generating functional is independent of the scale A we find 
the relation 



p * * Jp \ uifjp u «fp I 

+ / (|x- p <9aVp + \ ip-pQ^Xp +X- p Qa 1 ^p + iP-pQa 1 Xp) + (/a 



(2.7) 



where the anomalous dimension is defined to be 

V>p = |^p, ^p = |Vy (2.8) 

We can now combine Eqs. (2.6) and (2.7) to write, after some straightforward algebra, an equation for 
S). It will be satisfied if 



Pa(p)^-^Pa(p) 6S 



\84>p 5iP_ p Sip p 5tP_ p 

r ~i ( 2 - 9 ) 

+ 1 [Wp ( Pa {p) p ^ {p) ) +> - ip) Pa {p) ) W, 

Note that our claim is that if (2.9) holds for functional, as it stands, then a similar equation will hold for 
expectation values, which, in turn, ensures that all the Green's functions are invariant under the flow. We 
have just found, therefore, the sufficient condition we were looking for. 

We still have to take into account the effects produced by the rescalings needed after any Kadanoff type 
of change in order to complete a RG transformation. We have partially included them when we consider not 
bare fields but renormalized ones with some anomalous dimension rj. What is left over is just the canonical 
evolution of all quantities. To compute them in closed form the easiest way is to write (2.9) after having 
rescaled all variables with the appropriate powers of A to make them dimcnsionless, 

' SS . , 6S S . , SS \ 



S = 2(2*)'/*V)(#*#--7^ 

Jv \ Hp Sib_„ 0ip p 



P K H P Stp 8ip p Sip I 

(2.10) 



1-d + rj n2 K'(p 2 )\(- SS , SS\ f - „ d' SS , „ d' SS 



where the prime in d'/dp^ 1 means that the derivative does not act on the momentum conservation delta 
functions and thus only serves to count the powers of momenta of a given functional. 



Note that, once the above ERG equation is derived, it would be easy to derive a similar one for a 
model involving Yukawa couplings just by combining the present manipulations with that of, for instance, 
Refs. [6,38]. (The resemblance of (2.10) and a Polchinski type equation for scalar theories is pretty evident.) 

To fully specify the evolution of our theory under the RG flow, we have also to write down equations 
for the terms containing \v an( ^ X P an d the term with neither the sources nor the fields. Although we will 
not need them in our analysis, just for the sake of completeness we present the expressions obtained without 
further comments 4 : 

Q A (p) = P A (p) Q (p 2 ) , h = -J dt l <?~ 2 (P 2 ) X- A 1 (P 2 ) Xp , (2.11) 

with Q (p 2 ) a scalar function that evolves according to the equation Q (p 2 ) = |Q (p 2 ). 

As a final comment about our equation is that we present it on Euclidean space as it is customary in the 
field. For our purposes, however, there is nothing special about the Euclidean formulation, as finally what 
one obtains is just a set of relations among coupling constants. In fact, we have also derived the counterpart 
of Eq. (2.10) for Minkowski space. It is not so nice because of the extra presence of an imaginary unit coming 
from the functional derivatives of the Minkowskian "Boltzmann" factor e tS in the second term. Nevertheless, 
with this equation we have computed the /3-functions for a simplified action (one without operators with six 
fields) in much the same way we will explain later for Euclidean space: they are finite, real and consistent 
with the desired symmetries, as they should be. We have not proceeded further, but the parallelism between 
them and their Euclidean counterparts strongly supports the common lore that both should contain the 
same physical information and that the choice of space is much a matter of taste. Nevertheless, it would 
probably be nice to afford a complete calculation in Minkowski space. 

3. The action 

In this section we begin the discussion of an explicit example. We first define it through its symmetries, 
then justify how one can truncate its general action while still retaining non-trivial information and, finally, 
we give the prescriptions we have actually used to build it systematically. 

The sample model is that with N spin 1/2 two-dimensional Euclidean fields that obey the discrete 
symmetries of parity, charge conjugation and, to obtain reflection positive Green's functions, reflection 
hermiticity (see Rcf. [39] for a precise definition of them). We further impose the continuous symmetries of 
Euclidean invariance and the chiral symmetry U(N)r x U(N)l. 

For the definition to be consistent, one has to check that the above classical symmetries of the action 
will survive after quantization. That is, one has to ensure that the symmetries will be satisfied at any point 
of the flow if they are satisfied by the initial conditions. In our case this is verified nearly immediately by 
just looking at Eq. (2.10). The point is that the Kadanoff terms, which are the eventually dangerous ones, 
essentially take the form, in spinor and flavour indices, of the free kinetic term of the action. 

The next step is to choose an appropriate truncation. One would desire a kind of derivative expansion, 
at least because it is quite efficient when applied to bosonic theories [15,19]. However, the similarity with 

4 For a discussion with respect to the scalar case see again Ref. [38]. 
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the scalar case cannot be carried that far. The first important difference is that, unlike the scalar case, the 
zero momentum approximation (effective potential) is not feasible and the leading order is one with zero and 
one derivative terms. The reason is almost evident: Eq. (2.10) contains, due to the sum over polarizations, 
a j> factor in the Kadanoff terms, while a similar equation for bosons does not 5 . 

Another significant difference is that in the scalar case a general potential contains an infinite number 
of independent functional, whereas for finite N a general product of fermionic fields with fixed number of 
derivatives has in any case a finite number of terms due to the statistics. It is impossible to put twice the 
same Grassmann quantity at the same point. Thus, for the fermionic case the derivative expansion leads 
unavoidably to a polynomial approximation. This has practical consequences: the ERG equation becomes 
a large system of coupled non-linear ordinary differential equations instead of a small set of coupled partial 
differential equations, and the techniques to obtain numerical results arc different. Furthermore, the number 
of different structures for an arbitrary large value of N grows extremely fast as the order of the derivatives 
increases. In practice, it becomes practically intractable at order 3 unless the degree of the polynomial of 
the fields is also truncated. For this reason we work, up to a finite number of derivatives and also up to a 
finite number of fields. The remaining decision is to choose where to truncate. 

We require that a sine qua non property of a decent approximation is to allow a nontrivial anomalous 
dimension. Therefore, we will keep as many derivatives as needed to allow for a non-zero 77, within a 
reasonable number of fields (of the order of, say, twice the number of derivatives). With this criterion it is 
easy to realize that one derivative and four fields do not work: only the Gaussian fixed point is obtained, 
with classical critical exponents. Two derivatives seem in principle sufficient. However, once the /3-functions 
are obtained, it can be shown that the result r\ — is unavoidable, thus forcing us to work with terms up to 
three derivatives and six fields. 

The final preparatory step is to write down the action. To construct it systematically we list all 
symmetries and study the restrictions imposed by each of them. 

We will work with the momentum representation and, in order to simplify the notation, we will take the 
convention that any product of fields should be eventually integrated over the momentum carried by each 
field, with a common momentum conservation delta function. This would correspond to an integral over the 
whole space of a product of fields and their derivatives (of any order) at the same point. 

i) U(l). We begin with U(l), fermion number conservation. Its consequences are well known: the action 
must be built up of operators of the form 

sf 2 ^T(Pi)^ b (P2), p^ b = T(pi)^ b (P2), v? 2 ab ^T(Pi)i^ b (P2), (3.1) 

where we work in the momentum representation, a, b denote flavour indices and from now on the subindices 
of S, P, V label the fermion momenta. The Clifford algebra is defined by {7*,7 J } = 2<5 y with 7 S = — i^ 1 ^ 1 . 
Note that in two dimensions there are no other spinorial structures, since = i(? k ^ k . 

ii) Euclidean invariance . The Euclidean invariance is also easily taken into account: one has only to 
make sure that all Euclidean indices are properly contracted. 

iii) SU(N). The next one is the (vector) SU(N) group. If the fields transform under the fundamental 
representation, all possible scalar operators can be classified with the aid of Fierz reorderings. In fact, it is not 

5 See, for instance, Eq. (18) of Ref. [6]. 
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difficult to show by means of Fierz transformations that a general local operator in the trivial representation, 
built from products of fermionic fields, can be factored in terms of 

*12 = *12 7 ^12 = ^12 ! • / 12 = V 12 ■ \ 6 - Z ) 

Thus, the simplest manner to get rid of the internal group indices is to work with a basis written as products 
of scalar, pseudoscalar and vector operators (S, P, V^), transforming under the trivial representation, and 
powers of momenta. Therefore, the simplicity of two dimensions has come to help us again: a general 
functional can be written in terms of only three "building blocks" , and momenta. 

iv) SU R {N) x SU L (N). To enlarge SU(N) to SU R (N) x SU L (N) we realize that the chiral invariant 
operators are constructed from the combinations 

Vf2, £12634 — -P12-P34, S12-P34 — P12S34. (3.3) 

Therefore, if we restrict our attention to those kinds of terms, again with an arbitrary structure of momenta, 
the symmetry will be fulfilled. 

Note that the first type of operator in (3.3) carries a space index j and two fields, whereas the other 
two have no indices and four fields. From this it can be immediately inferred that with an even number of 
derivatives one can only have operators with An fields (n integer), whereas with an odd number of derivatives 
the allowed operators contain An + 2 fields. The reason is that all indices must be contracted, either by the 
Kronecker delta tf*- 7 ' or by the complete antisymmetric tensor in two dimensions e*- 7 ', which implies that an 
odd number of derivatives needs an odd number of operators of the type V{ 2 . 

v) Parity. It only remains to impose discrete symmetries. Parity is easy: for a Euclidean invariant 
operator, products of S, and momenta are parity-conserving. The only problem is when we have the 
pseudo-scalar operator P. What we have to do is just follow the standard rule: a term with an odd number 
of P's must contain a Levi-Civita symbol e y also; a term with an even number of P's must not. 

vi) Charge conjugation. To impose charge conjugation and reflection hermiticity proves to be the most 
involved task. This is because both operations exchange fermions and antifermions, and thus they change, in 
general, the momentum structure. Explicitly, under charge conjugation our elementary operators transform 
as 

Si2->S 2 i, Pi2^~P2i, (3-4) 

To take into account this symmetry at the level of the basis, the most effective manner is to consider all 
momenta written in combinations like (pi ±p2) J , where p\ is the momentum of an antifermion and pi the 
momentum of the fermion of the same bilinear. In this way it is easy to distinguish between C-conserving 
and C-violating operators, and to construct both sets. 

vii) Reflection hermiticity. The last one is reflection hermiticity It is defined, in principle, in coordi- 
nate space [39] and under such transformation, our "elementary operators" behave just as in (3.4). What 
is new is that when transformed, one must change the coefficient of the operator by its complex conjugate. 
Therefore, once we restrict ourselves to C-conserving terms, this additional symmetry restricts the coeffi- 
cients of those terms to be real. The only subtlety is that, as it is defined, the fields do not become complex 
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conjugated and neither do their derivatives 6 . And if one remembers that a derivative in coordinate space 
amounts to a factor — ip this indicates that an extra power of i should be added for each power of momentum. 

viii) Further degeneracies. Finally, the freedom of integrating by parts (each operator has a delta 
function of momenta conservation) relates different functionals, and, ultimately, reduce the number of in- 
dependent ones. The best way of implementing these final constraints is to find out a criterion in order to 
write down every operator in a "standard" way. We will explain ours in Appendix A, where we will also 
write down the complete action, consisting of a basis of 107 functionals. 

One of them is rather peculiar. It is 

iV( 2 (p 1 -p2)Hpi-P2) 2 - (3.5) 

One may be worried about it because it would lead to a propagator with an additional pole besides the 
physical one on the particle mass-shell, thus entering in conflict with unitarity. This is, however, not 
important at all, because the above kind of reasoning implies that one assumes a well-defined perturbative 
expansion, and this is not the case (we have irrelevant operators that make any perturbative expansion 
around the Gaussian fixed point completely ill-defined). One should think that the theory is such that it 
manages to have a well-defined complete two-point function free from unphysical singularities. A completely 
different point is that, besides the above discussion, when one computes the /3-functions of the theory one 
realizes that this operator, at least up to the order we are considering, does not contribute to any other. 
Therefore its evolution will affect absolutely no conclusion we obtain without it. For this reason, we do not 
include it in the action. We should remark, however, that for all we do it is as if this term were already 
there, although for the sake of brevity we will not write it down any more. 

4. Computing the /3-functions 

Once we have constructed the initial action we want to work with, the next step is to substitute it into 
the ERG equation (2.10) and to compute the /^-functions of our model within the given approximation. 

In principle this is just a purely algebraic exercise. Nevertheless, it turns out that from a practical 
point of view it is an almost forbidding task, if done by hand. During intermediate steps of the calculation 
one has to handle thousands of terms and it is too easy to make errors. For instance, when computing 
j^-Pa (p) J^ S , the functional differentiation gives 302 terms, and one has, roughly speaking, to square 
them and multiply the result by the inverse propagator. Then, one has to compute the appropriate products 
of gamma matrices, expand all the terms and, finally, perform the integration by parts to reach our chosen 
basis. The number of operators considerably increases in these last processes. Thus, it is mandatory to use 
a symbolic manipulation computer program to perform the functional differentiation, do the algebra and 
integrate by parts. Because of this, our computation was done with the help of Mathematica. 

To calculate the flow equations, we use an extended action, greater than that discussed so far, in order 
to have some extra check of our equations. That is, we consider an action expanded in a basis that consists 

6 We remind the reader that, in order to turn properly from Minkowski to Euclidean space, one has to redefine 
the symmetries of the problem, specially those which involve complex conjugation. Our definitions coincide with 
those of Ref. [39]. 
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of terms with two fcrmions with one and three derivatives, four fermions with zero and two derivatives and 
six fermions with one and three derivatives, but without imposing any symmetry other than vector U(N), 
parity and Euclidean invariance (that is, we impose neither reflection hermiticity, nor charge conjugation, 
nor chirality). We then project the space generated by this basis into the invariant subspace under the 
required symmetries and its direct complement. The required flow equations are obtained after the first 
projection, while the complementary subspace provide us a consistency check of the calculation. They define 
a set of null equations that have to be satisfied along the renormalization flow: after projecting to an initial 
symmetric action, any non-zero contribution of a non-symmetric term will indicate an anomaly, which we 
have argued are non-existing. We leave to Appendix B the complete set of /3-functions. 

Finally let us justify the inclusion of operators with three derivatives into the action. As we advanced 
above, it is motivated by the fact that those terms arc necessary in order to get a non-vanishing critical 
anomalous dimension. The argument is as follows. The anomalous dimension is related to the fact that we 
are free to fix the normalization of one term of the action, by choosing an appropriate normalization of our 
fields. If, as is customary, we keep fixed the coefficient of the so-called kinetic term, then its /3-function is 
substituted by an equation for n which, in practice, is calculated in a similar fashion. We have to study, 
then, which type of Kadanoff transformations contribute to ip_ p fiip p . The term j^-ifi JL S cannot, unless 
there were a mass operator, which is forbidden by chiral symmetry. There are, however, some contributions 
coming from S^ 2 \ due to j^-ifi JL S . We will find, hence, the anomalous dimension as a linear combination 
of couplings of S^ 4 ' 2 ) and, consequently, if these couplings vanish at the fixed point then r\ = is unavoidable. 
(We define <S( a ' b ) as the part of the action that contains a fields and b derivatives). If one now studies S^ 4 - 2 \ 
it is not difficult to convince oneself that its only contributions must come from S^ 6 ' 3 ' , apart from canonical 
rescalings. The implications are now immediate: if S^ 6,3 ) did not exist, then the whole action <S>( 4 ' 2 ) would 
evolve canonically, thus it would vanish at the fixed point and we would obtain a vanishing anomalous 
dimension. 

5. Fixed points, critical exponents 

5.1. Generalities 

The next step is to find the fixed point solutions, that is, the sets of coupling constants that make all 
the /3-functions vanish. These will indicate the points to which the RG tends to, thus providing us with the 
first indication of what the phase diagram of the system looks like. 

The condition S = is equivalent to a system of 106 non-linear algebraic equations. To simplify it we 
note that all the coupling constants of operators with six fields must enter linearly, because the only source of 
non-linearity of Eq. (2.10) is its first term on the r.h.s., and it can give contributions neither from S^ 6 ' 1 ^ nor 
from iS^ 6 ' 3 ), within our approximation. Therefore, we can reduce the system to a set of only 13 non-linear 
equations, 

= 277 5l + 8glajN/(-2 + 3r?) + {8glf35N + 8/3 7 [<?i(-4r 2 + si - s 2 - 3s 3 - 3s 4 ) 

+ 2#2 (mi - m 2 - m 3 ) + 40im 3 iV]}/(-4 + 3jj), 
= 2r7.g 2 + 83?a 7 /(-2 + 3t?) + {8. 9 2 /3<5 + 8/3 7 [2 5l (mi - m 2 + m 3 - s 2 ) 
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+ 5 2 (-4r 2 + si -8 2 - 3s 3 - 3s 4 ) + 2g 2 s 2 N]} / (-1 + 3r?), 
= 2(-l + 77)7711 + {16.gia(5(-2.g2 + 9\N) + 4a~/[gi(mi ~m 2 ~m 3 - 2r 1 - 6r 3 - 4s 2 - 3s 3 + s 4 - i) 

- .g 2 (3mi + m 2 + 5m 3 + 2t) + 2 gi N(2 mi + 3m 3 + r 2 + 2s 3 + s 4 + t)}}/(-4 + 3r>), 
= 2(-l + r/)m 2 + {16ry lff2 a5 + 2a 7 [ 5 i(2ri + 2r 2 + 6r 3 - 3s a + 5s 2 + s 3 + 5s 4 ) 

+ 2 32 (mi +m 2 + 3m 3 - i)]}/(-4 + 3r?), 
= 2(-l + 7?)m 3 + {16gi.g 2 CH5 + 2a7[ 3 i(2mi - 2m 2 - 2m 3 + 2r 4 - 2r 2 + 6r 3 - si - s 2 - 7s 3 - 3s 4 - 2i) 

+ 2 32 (2m 3 + 2m! - t) + A gi N{m 3 + r 2 + 2s 3 + s 4 )]}/(-4 + 3r?), 
= 2(-l + r/)n + {8g 2 a5(2g 1 - g 2 N) + Aa-f[gi(mi - m 2 + m 3 + 4r 3 - 2s 2 ) 

+ fl2(3mi + m 2 + m 3 - 2r 2 + s a - 3s 2 - 4s 3 - 4s 4 -t) + 2g 2 N(r 2 - 2r 3 + s 2 + 2s 3 + s 4 )]}/(-4 + 3r>), 
= 2(-l + V )r 2 + {A(g\ + Agl)a5 + 4a 7 [ 5l (2m 3 - t) + .92(277 + 6r 3 - Sl + 3s 2 )]}/(-4 + 3r?), 
= 2(-l + r?)r 3 + {12^a5 + 4a7[gi(3mi +m 2 + 5m 3 - t) 

+ 32(^1 -m 2 -m 3 - si - s 2 - 3s 3 + s 4 - t) + 2g 2 N(r 2 + 2s 3 + s 4 )]}/(-4 + 37/), 
= 2(-l + r/)si + {16g 2 a8(2gi - g 2 N) + 16a~/{gi(2r 3 + s 2 ) + g 2 (mi + m 2 + m 3 ) 

- ff2 7V(2r 3 + s 2 )]}/(-4 + 3r ? ), 
= 2(-l + 77)s 2 - {8#?a(5 + 8.9107(7711 + m 2 + m 3 )}/(-4 + 3rj), 

= 2(-l + 7/)s 3 + {Ag\a8 + 4a"f[ gi (2 mi + 2m 2 + t)+ g 2 {-2r x + 2r 3 + Si + s 2 )]}/(-4 + 3r?), 
= 2(-l + 77)54 + {4c,2 a( S + 4a7[ 5 i(2m 3 - t) + g 2 {2n - 2r 3 - 8l - * 2 )]}/(-4 + 3t?), 

= 2(-l + 7?)t + 8a7[3i(-2n + 2r 3 + s a + s 2 ) + .92(7711 + 3m 2 - m 3 ) + .9i7V(-2m 2 + t)]/(-4 + 3r/), (5.1) 
where 77 is the anomalous dimension that turns out to be 

rj = 4a[-m 4 + m 2 + m 3 + si + s 2 + s 3 + s 4 + t - 2N{r 2 + 2s 3 + s 4 )]. (5.2) 
N is the number of flavours and a, (3, 7, S are scheme dependent parameters defined as 

a =J^J p K '(P 2 )> P = J^y j/ K> (P 2 ) ' 7 = ^(0), f = *"(0). (5.3) 

The appearance of the above quantities just reflects the freedom in choosing a renormalization scheme. 
Furthermore, although the /3-functions depend on four parameters, we will see that the fixed point solution 
will depend only on two combinations of them. This is just the pattern that occurs in a scalar theory within 
a similar truncation [19]. 

The system cannot be solved analytically, unless we perform further simplifications, like keeping only 
the dominant term in an asymptotic expansion at N —> 00. On the other hand, one can, of course, simply 
try to study its solution numerically. We will present both approaches in turn. 

After the fixed points are identified, the behaviour of the theory near each of them is controlled by the 
critical exponents. One of them is fixed once we solve our set of equations: it is the anomalous dimension 
at the fixed point value. The rest are found by linearizing the RG transformations near the chosen fixed 
point. That is, if gi is a generic coupling constant, then its variation in the vicinity of a fixed point go 
is approximated by Sgi = gi = Rij\ ga 5gj, where Rij is the matrix The eigenvalues of Rij\ go can be 
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identified with critical exponents. They can be thought of as the anomalous dimension of the operators 
which drive the theory away from the fixed point. 

We can now no longer work with the reduced system of 13 couplings, but the full 105 x 105 matrix is 
needed as we allow deviations from the fixed point values of the six-fermion couplings. Around the Gaussian 
fixed point, for instance, the four-fermion and two-derivative operators have the same degree of "irrelevance" 
as the six-fermion and one-derivative ones. 

Finally, let us turn again to the issue of scheme dependence. We have just said that, in general, the 
precise values of the coupling constants at a fixed point are scheme dependent, thus reflecting that they are 
not universal quantities. Critical exponents, on the other hand, are universal, hence they should be scheme 
independent. Nonetheless, due to the truncation, scheme dependencies will inevitably appear. What we will 
do is, as usual, to try to find a suitable scheme where the dependence will not be that important. To this 
end, we will apply to the various solutions a principle of minimal sensitivity to discriminate among different 
results [40,19]. 

5.2. N -> oo 

We are now going to set up a large N expansion for our model, with which analytic results can be 
obtained. Later on we will see that when we study the general case by suitable numerical approximations, 
we will recover our present results as a first term of the asymptotic series around N — > oo. 

To define properly our approximation, we substitute each coupling constant gi by N Zi gi and study the 
limit N — ► oo keeping gi fixed. In principle, Zi can be any real number, but for simplicity we only consider 
integer values. We then find the set {z^ that makes all /?- functions finite and, if possible, non-trivial. 

With the above requirements, there are essentially two different manners to define the 1/N expansion, 
which lead to different results. We label them by I and II, and discuss each in turn. 

The Type I solution is obtained by considering Zi = — 1, where i runs over every of the couplings that 
enter in Eq. (5.1). With this definition, the anomalous dimension vanishes at leading order in 1/N and the 
system (5.1) becomes 

= -4a75i - 2/3<Sg 2 - SI3^gim 3 = -Ag 2 s 2 l3l, 
= -2mi - AaSgl - 2ajgi{2m 1 + 3m 3 + r 2 + 2s 3 + s 4 + t) = 
= -2n + 2a5g\ + 2ajg 2 {-r 2 + 2r 3 - s 2 - 2s 3 - s 4 ) = -2r 2 
= -2si + Aa5gl + 4a~/g 2 (2r 3 + s 2 ) = -2s 2 = -2s 3 = -2s 4 , 
= -It + 2ajg 1 (2m 2 - t). 

Its solution is 

g 1 = -l/(aj), rn 2 = 0, m 3 = (5/(4a7 2 ) + 1/(2/37), n=aSg 2 , r 2 = r 3 = 0, 
s 1 =2a8gl, s 2 = s 3 = s 4 = 0, t = 5<5/(4a 7 2 ) - 6/(4/3 7 ) - m x . (5.5) 

We now choose = —2 for all the six fermions coupling constants. This is not the only solution since 
there are other rescalings consistent with the reduced system (5.4). For example, one can assign the value 
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-2m 2 = -2m 3 - 2ajg 1 (m 3 + r 2 + 2s 3 + s 4 ), 
= -2r 3 - 2ajg 2 (r 2 + 2s 3 + s 4 ), 



(5.4) 



— 1 to some of the Zi's and —2 to the other ones, but it turns out that the results below do not depend on 
that. 

The characteristic polynomial P(X) associated to the matrix of linear deviations is exactly computable, 
P(X) =A 2 (A + 2) 12 (A + 4) 83 (A + 6) (A 2 + 6A - 8) 

(5.6) 

x (-A 5 - 12A 4 + (8w - 44)A 3 + (64w - 16)A 2 + (32w + 64)A - (128™ + 256)) , 

where w — (35/(aj) and z = 5/j 2 . We can read the critical exponents from P(A). There are 100 scheme 
independent eigenvalues, most of them coinciding with the Gaussian values 0, —2, —4 and —6. The non-trivial 
ones are —3 + vTf = 1.1231... and —3 — \/l7 = —7.1231..., and the roots of the polynomial 

Q(A) = -A 5 - 12A 4 + (8w - 44)A 3 + (64w - 16)A 2 + (32w + 64)A - (128w + 256), (5.7) 

which are all w-dependent. If w < 0, that corresponds, for instance, to the exponential cut-off function 
K\(p ) = e s a i ! the more relevant critical exponent is Ai = 1.1231... 

Note that the fixed point solution depends freely on g 2 and m\. This is the expected result for the 
chiral Gross-Neveu model because the £7(1) Thirring like excitations (which in our action are controlled by 
g 2 ) decouple from the rest and this subsystem is conformally invariant (i.e. it is at fixed point) for any value 
of g 2 . For the SU(N) part there exists a discrete set of fixed points, the one of Dashen and Frishman being 
one of them. This critical point is reached when the constant g\ is of order l/N, as in our case. So we can 
make a correspondence between our solution and that of Ref. [35]. Nevertheless, the values of the anomalous 
dimension in both cases do not match. For the cited fixed point it is non- vanishing at leading order in l/N, 
and not zero as we have found. This discrepancy with the exact result of Dashen and Frishman could be 
caused by our truncation. We cannot reject, however, the possibility of having found a different fixed point 
as it has already occured previously [36] . 

For the Type II solution it is useful to define the new variable vn\ — mi — m 2 + m 3 instead of m\. Then 
it corresponds to the following rescaling of couplings 

gi^gi/N, g 2 ^>g 2 /N, m / 1 —>m' 1 , m 2 —>m 2 /N, m 3 ^m 3 /N, n^n/N, r 2 -> r 2 /N, r 3 ^r 3 /N, 
si -» si/TV", s 2 -» s 2 /N, s 3 -> 83/N, s 4 -> s 4 /N, t -» t. (5.8) 

The solution takes the form 

2a-/?7 -2a + /3 7 

.91 =ff „ 3 > 3i= P- 



48a 3 ' " ' 24a 3 

-8a + /? 7 8a - /3 7 24a 2 (35 + 32a 3 7 - 18a/3 2 <5 7 - 18a 2 /3 7 2 + 3/3 3 c5 7 2 + 4a/3 2 7 J 



mi " 144a 2 ' m2 ~ 72a 2 ' ™ 3 " 576a 3 7 (-4a + (3j) 

_ 43 (3 2 5 (35 43/3 7 -8a + /3 7 _ 8a - /3 7 

ri ^36^ + 12^3 ~ 6a^ ~ 288a^' r2 ~ 288a 2 ' ^ ~ 36a 2 ' 

20 (3 2 5 (35 5/3 7 8a - /3 7 -8a + /3 7 -8a + /3 7 

s l = « 1" 7 q _ 7T^ "To — 9 ' s 2 = : , , 9 , S3 = — — , S 4 = 



'9a 6a 3 3a 2 7 18a 2 ' " 144a 2 ' J " 288a 2 ' 4 " 288a 2 
8a - /3 7 
144a 2 



which has a non-zero anomalous dimension, 
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z = 0.1 


z = 0.5 


z = 1.0 


z = 2.0 




2.258 


2.239 


2.217 


2.175 


w* 


0.122 


0.616 


1.250 


2.610 


V 


1.130 


1.128 


1.125 


1.116 



Table I: Local minimum of Ai, the most relevant critical exponent, corresponding to 
the solution for N — > oo labelled as Type II in the text, for different values of z. n is 
the anomalous dimension at that point and w* is the value of the parameter w at which 
the minimum is reached. 

The set of the remaining z^s is unique and composed of the numbers —1 and —2. Unfortunately, 
unlike the previous case we could not find the exact analytical expression of the characteristic polynomial. 
However by computing numerically the eigenvalues for different values of z and w we could guess some exact 
results. None of the critical exponents coincide with their canonical counterparts. Moreover, most of them 
are functions of the combination ^. Thus there are 82 eigenvalues A = — |j, 8 eigenvalues equal to |(1 — f^) 
and 4 of the form 2 — The remaining ones are not functions of the ratio w/z only (and even a few have 
a non- vanishing imaginary part, which is not unusual in approximations based on truncations). We have to 
study numerically the most relevant critical exponent, which belongs to the class with no simple dependence 
in w and z, for different scheme paramctrizations. As it happens in the scalar case, for any value of z, this 
exponent always has a minimum at some scheme parametrized by w = w* . This behaviour induced us to 
use the minimal sensitivity criterion to fix the parameter w to its critical value w*. Unfortunately due to 
the monotonic dependence of the solution on the parameter z in the range analyzed, we were unable to set 
it with a similar prescription. In Table I we show some values of A* = \i(w*) and the anomalous dimension 
for various z's. 

5.3. Finite N 

For a finite number of flavours analytical results for the fixed point couplings cannot be found. So one 
has to proceed numerically to search for the zeroes of the /^-functions. Moreover, the number of different 
solutions of a system of coupled non-linear equations is not known a priori and the common routines for root- 
finding (such as the FindRoot command of Mathematica) do not guarantee that all the zeroes are reached. 
A more serious inconvenience is to decide if a given zero corresponds to a real fixed point solution or if it is 
a spurious root resulting from our truncation. 

The first problem can be reasonably reduced after some experience is acquired. In fact, we can know 
by intuition which is a reasonable range of values for the couplings and inspect this region exhaustively. Of 
course this is not easy for a system of thirteen equations, but we can gain some confidence in the results if 
we examine minutely the adequate region. 

The second problem, however, is much more complicated. In principle we do not have any criteria to 
decide if a root of the /3-functions system corresponds to a genuine fixed point solution or if it is a fictitious 
artifact of our approximations. This problem, which already appeared in the bosonic case too, is perhaps 
the Achilles' heel of the approximations based on truncations [17,18]. We present the class of solutions of 
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which we are more confident. These are mainly the ones which asymptotically match with some solution 
clearly identified in the framework of the large TV expansion. 

Let us, to begin with, select a particular scheme and find the solution for different values of TV: w = —2 
and z = 0.5, corresponding to an exponential regulating function (K(x) = e~ x ). We analyze the dependence 
on these parameters later on. 

One solution is found which behaves asymptotically as the type I one in the TV — > oo limit. TV 77 increases 
with TV and tends to 4.87..., while the most relevant critical exponent Ai decreases with TV asymptotically to 
the value 1.1231..., in agreement with the 1/TV expansion. For the second eigenvalue we find complex figures 
that we attribute to our approximations. Another piece of bad news is that, unlike the TV — > oo case, the 
solution for finite though big enough TV, is isolated, while, as we mentioned before, the fixed point solutions 
for Thirring like models are continuous in the U(l) sector. Again we blame this confusing result on the 
truncation. We present in Fig. 1 the curves for Nrj and Ai as a function of TV. 

More interesting is perhaps the study of the dependence of the solution on the scheme. We have noticed 
that z enters the equations only through the anomalous dimension as a global factor. For this reason, the 
dependence of the fixed points solutions in z is quite simple: it is almost linear in rj. Therefore it is more 
attractive to investigate its behaviour under changes on the parameter w for fixed TV and z. The motivation 
of this analysis is the search, as in the scalar case, of some non-linear w-dependence in such a way that we 
can invoke a principle of minimum sensitivity to fix the value of this parameter and eliminate one fictitious 
dependence. To this end, we fix the value of z to z — 0.5 and TV to TV = 1000. The curve rj vs. w is 
monotonic and decreases with w, while the first eigenvalue Ai reaches its minimum value Ai = 1.12511 at 
w = —45, which increases as we lower TV: it is equal to 1.1273 for TV = 500 (it is reached at w = —23), 1.146 
for TV = 200 (at w = -10), 1.1519 for TV = 100 (at w = -8), 1.695 for TV = 10 (at w = -2.4) and finally, 
2.560 for TV = 3 (at w — —0.5). We show two of these curves in Fig. 2. 

For the fixed point that matches the Type II solution as TV — ► oo we found a curious behaviour. For 
iV moderately large, (say TV = 1000), the numerical solution is in good agreement with the 1/TV analytical 
result (for example, the value of r\ for z — 0.5 and w — — 2 is 1.99, compared with the exact -q = 2 for 
TV — > oo). As we lower TV the values of the anomalous dimension rj and the most relevant eigenvalue Ai 
decrease. But, unexpectedly, the solution disappears at TV = 142 (actually at TV — 142.8 if we let TV take 
non-integer values). A closer analysis of the space of solutions shows us that at this value of TV the branch 
of solutions compatible with the type II 1/TV expansion merges with another family of fixed points. This 
last branch has finite asymptotic limits for rj and Ai as TV — > oo. However some couplings do not behave 
as a power of TV in this limit and, therefore, it cannot be associated with a 1/TV fixed point in the sense 
stated previously. At the bifurcation point rj = 1.88 and Ai = 5.80. We show in Fig. 3 the curves rj(N) and 
Ai(TV). This peculiar behaviour of the type II solution suggests that it cannot be identified with the Dashen 
and Frishman fixed point, which exists for any value of TV. Even though this behaviour could be another 
consequence of the truncation, it is hard to justify it because for low N only operators with few spinor fields 
are allowed by the Pauli principle, and thus we expect our six-fcrmions truncation to be accurate. We have 
not been able to solve this puzzle. 

It is easy to find many other solutions, especially for low N. For some of them there exists a minimum, 
either for r/ or for Ai but in other cases both curves are monotonic in w. They have also different behaviours 
as N — > oo. We show one example in Fig. 4. 
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54- N = 1 

Finally we will consider the special case N = 1. For this particular value of N not all the operators 
presented in Appendix A are independent. As we mentioned in Sect. 3 the effect of the Fierz transformations 
is to relate covariant U(N) local operators (like ip a (p\) ip b (JP2)) to scalar ones (like ip a (pi)ip a {P2))- So in 
the N = 1 case the Fierz transformations uncover relations between the S, P and the operators. For 
example, for operators without derivatives we have the identities 

Si 2 S 34 = -P 12 P 34 = —V^Vl. (5.11) 

They establish relations between the coupling constants that permit us to reduce considerably the 
system. For the set of couplings of the four fermions operators, the independent ones are 

9 = .91 - 92, u 1 =m 1 -m 2 + m 3 -r 1 +r 2 -r 3 + (si - s 2 + s 3 + s 4 )/2, 
u 2 = 2mi + 2m 2 - 2m 3 + 2r\ + 2r 2 - 2r 3 , ii 3 = 4m 2 + 2n - 2r 2 - 2r 3 - 2s 3 - 2s 4 , 
U4 = —81 - s 2 + s 3 + 84, U5 = -2s 3 + 2s4 + 2t. (5-12) 

Eq. (5.1) is now a 7-equation system that looks like 

= 2g(4rj - 3rj 2 + 4u 2 w - Aii 3 w + 4u 4 w)/(4 - 3??) 

= 2(8g 2 + 4ui + 8giii + 2gu 3 + 2.gu 4 - 7ui7? + 3uirj 2 ) /(-4 + 3rf) 

= 2(8g 2 + 8gu! + 4u 2 + 8gu A + 4gu 5 - 7u 2 i] + 3m 2 ?7 2 )/(-4 + 3j?) 

= 2(-24g 2 - 24gu! + 4gti 2 + Au 3 - I2gu 3 + Agu b - 7u 3 ri + 3u 3 r] 2 )/(-4 + 3j?) 

= 2(8g 2 + 8gui + 4gu 3 + 4u 4 + 4gu 4 - 7u 4 r] + 3sry 2 )/(-4 + 3rf) 

= 2(-4gu 2 - 4gu 3 - 8~gu A + 4u 5 - 7u 5 rj + 3m 5 ?7 2 )/(-4 + 3ry) 

77 = — 2u 2 z + 2u 3 z — AH4Z + 2u 5 z. (5.13) 

It is linear in u\, u 2 , u 3 , 114 and U5, so we can solve it for these variables ending with a two-equation system 
for g and 77. After a bit of algebra and discarding the trivial solution we finally get a unique equation for 77, 

= -120w 2 z + 288wz 2 + r)(13w 2 - I32wz + 210w 2 z + 288z 2 - 720wz 2 ) 

(5.14) 

+ rj 2 (99wz - 90w 2 z - 432z 2 + 594tjjz 2 ) + r/ 3 (162z 2 - 162w;z 2 ). 

As in the previous analysis we have to choose some particular scheme, i.e. fix w and z, and solve the 
equation numerically. Unfortunately, a simple inspection of the equation reveals bad news. The system is 
not undetermined and there is no room for a free ^-dependence of the fixed point solution as it is true in 
the Thirring model. This property is satisfied in the previous order approximation (terms with less than 
three derivatives) where r] vanishes identically. The reason this property is lost in the three derivatives 
approximation is unclear for us. For a more detailed analysis it is necessary to go to the next order to see 
whether this property is restored. 

We solved Eq. (5.14) numerically for different values of w and z. As in previous examples the fixed 
point solutions are almost linear in the parameter z so it is more interesting to study the behaviour of the 
solution as a function of w. However, in the range of values studied, we did not find any non-monotonoic 
behaviour either in the critical couplings or in the anomalous dimension. We present some of the results in 
Table II. 
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w = -0.1 


w = -0.5 


w = -1.0 


w = -2.0 


z = 0.1 


1.763 


3.691 


6.316 


11.747 


z = 0.5 


1.418 


1.790 


1.418 


3.388 


z = 1.0 


1.376 


1.559 


1.811 


2.349 


z = 2.0 


1.354 


1.445 


1.569 


1.834 



Table II: Values of r\ for N = 1 and different scheme parameters z and w. 

6. Summary and discussions 

In this article, we analyze the application of the ERG method to fermionic theories. An ERG equation for 
Grassmann variables is derived and the critical properties of the chiral Gross-Neveu model in two dimensions 
are studied with it. 

To solve the ERG equation, a non-linear functional differential equation, we perform a double truncation, 
in the number of derivatives (derivative expansion) and in the number of fields (polynomial approximation) . 
Unfortunately, these approximations produce similar problems that already appear in the scalar case within 
analogous truncations: spurious solutions and unphysical scheme dependencies. The latter, which is a 
common feature of almost any approximation in QFT, can be partially disentangled by invoking a minimum 
sensitivity criterion: for a given observable we choose the scheme that gives the most "stable" result. The 
emergence of spurious solutions is a more serious problem. In principle we do not have any strong argument 
to accept or reject a solution, except for those which lead to absurd results. 

Note that for the bosonic case, within the derivative expansion, we can either expand the action as a 
polynomial in the fields, leading to a system of coupled non- linear equations or not make any further ap- 
proximation and consider the potentials as arbitrary functions (not necessarily real analytical), that requires 
the study of partial differential equations. While the first approach produces lots of spurious solutions of the 
fixed point equations [17,18], the former has shown to produce the correct ones [15,19]. For the fermionic 
case, however, the situation is quite different. For finite N, within the derivative expansion, a truncation in 
the number of fields is not an approximation for local Lagrangians, but the definition of a function in terms of 
Grassmann variables. So, for fermions, the polynomial approximation should not produce fictitious solutions 
if we are constraining the number of derivatives. In accordance with this, for N = 1, which is the only case 
where we actually have a pure derivative expansion without any further truncation, no spurious solution 
appears. (There are only three solutions, besides the trivial one, and two of them have complex coupling 
constants, thus being rejected at once). It would be interesting to perform a pure derivative expansion for, 
say, N — 2 and check if the above pattern holds. 

The first analysis we do of the fixed point structure of the chiral Gross-Neveu model is in the large 
N expansion of the /3-functions. We find that it can be defined in two different ways, with remarkably 
different results. The first one leads to a continuous family of fixed-points which reminds that of Dashen 
and Frishman: the solution is free in the direction associated to the abelian degrees of freedom and fixed of 
order 1/N in the direction of the SU(N) ones. However it presents an important difference: the anomalous 
dimension vanishes at leading order in contrast to the order 1 value of the Dashen and Frishman solution. 
We attribute this difference to the truncation. The inclusion of more terms should clarify this point. 
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The other type of solution is much more involved. Its anomalous dimension is non-zero but the Thirring 
like excitations do not appear. Moreover, unlike the preceding case its dominant eigenvalue of the linearized 
RG transformation depends on the scheme. Another astonishing result is the structure of the remaining 
eigenvalues. One would expect that the most irrelevant ones would not be too different from the canonical 
ones due to our truncation, which is not the case. Furthermore, the stability of the solution for finite N 
seems to indicate that it is not an artifact of the truncation, but a true fixed point. 

To go further we proceed numerically. We can clearly identify a solution that asymptotically matches 
the first of the above ones and trace it to very low values of N. It presents two important drawbacks. The 
first one is that it is isolated, unlike the strict limit N — > oo, where a one-parameter space of solutions 
appears. The other one is again a remaining dependence on a parameter which labels different schemes, 
although an accurate analysis of the most relevant critical exponent exhibits minimum sensitivity to some 
schemes. A search of the behaviour of the critical exponents as a function of N clearly shows that the value 
of Ai decreases with N while N-q increases. 

We can also find another set of solutions, for N > 142 that matches the second type of the large N 
ones. At N = 142.8 it merges with another family of fixed points, with divergent /3-functions when N — > oo, 
although its critical exponents seem to be finite in that limit. This odd behaviour ruled out an identification 
of this solution with that of Dashen and Frishman. The lack of information in the literature (anomalous 
dimensions, critical exponents) about the extra fixed points of the model does not allow us to recognize our 
solution as any of them. 

Finally, due to its peculiarities, we separately analyze the N = 1 case. The results are, however, 
discouraging. On one hand, at first order in the derivative expansion, there appears a solution with a free 
parameter, which labels the Thirring like excitations, but it gets spoiled when higher orders are considered. 
On the other hand, we find a severe two-parameter scheme dependence, which made unreliable any conclusion. 

Let us remark that the insufficient non-perturbative studies (lattice computation, etc.) of the chiral 
Gross-Neveu model in d = 2 impede to discriminate definitively in favour of our results. However we are 
very confident of some of our findings: as we argue above, within the 1/N expansion the first fixed point 
solution is an excellent candidate for the Dashen- Frishman fixed point, whereas the other one presents 
evidences to be a new fixed point, with quite intricate properties, not discussed previously in the literature. 
Moreover, both solutions have a smooth behaviour for finite N. 

We want to pause here to make a comment about redundant operators, that is, those operators that 
do not affect correlation functions [11]. Examples are functionals proportional to the equations of motion. 
They are known to be present in general and, specifically, one should expect the existence of a redundant 
operator that reflects the freedom of changing the normalization of our action (e.g. modifying the coefficient 
of the kinetic term). Its associate eigenvalue is scheme dependent and of no physical relevance. However, 
the ERG equation we consider presents an invariance under the normalization of the fields [14] which is 
not preserved by our approximation (not even in a pure derivative expansion). The signal of its restoration 
must be an appropriate redundant operator with vanishing eigenvalue. Its appearance should, fortunately, 
fix the anomalous dimension close to its correct value. It should also be mentioned that the presence of this 
operator has been used several times in the literature to discriminate among different schemes and find the 
appropriate anomalous dimension [41,16]. We expect that similar techniques should apply for fermions also, 
although we have not gone through the details. More work has to be done in this direction. 



18 



As a summary, we may say that, in general, our results seem somewhat discouraging, specially those 
for low N. Let us note, however, that the two dimensional world is rather peculiar, due to the importance 
of quantum effects, which generally produce large anomalous dimensions. Thus, technical simplicity turns 
into increasing complexity while the dimension is lowered. Nevertheless, we have gone much further than 
similar computations for bosonic theories, where in d = 2 it seems that the method completely breaks down 
[17]. Another interesting feature is the seemingly good results for the large N limit. At the computational 
level this is related to the fact that at the leading order of the 1/N expansion the system (5.1) of the flow 
equations simplifies dramatically and no room is left for spurious solutions. Actually, improvement of results 
in the large N limit appears to be a general feature of the ERG approach (see, for instance, Rcf. [10]). 
Therefore, the credibility of the method can only be clearly decided after extensions to other dimensions 
and, possibly, the inclusion of higher order terms. 

A last comment is dedicated to further work. As we have just mentioned, the formalism should be 
extended to higher dimensions. Equation (2.10) is prepared for that. What has to be done is to choose 
an appropriate action and perform a similar calculation. The number of spinor structures will be increased 
and, therefore, one will have to handle more terms in the action. However, we guess that, due to the 
greater complexity, two derivatives may be sufficient to obtain interesting results, or at least, according to 
the standard rule that quantum effects become less important when the dimension is increased, we hope 
that, within the same approximation, the results will be more transparent. 
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Appendix A. Conventions and notation 

In this appendix we summarize our notational conventions. 

We are working always on Euclidean momentum space. A standard term in the action should be 
/ [ ■ ■ [ f TT^ S pi+p[+---+Pn+P^ 1 i r ^ •••^r„Vp"x(polyn. in momenta) x (flavour matrices), 

Jp! Jp[ Jp n Jp' n \ Z7T ) 1 

(A.l) 

where 

j p ^j dd P> 6 Pr+P' 1 + -+P-+P' n = W*" (l> +P^j > ( A -2) 

the fermions are denoted indicating its momentum label as a subindcx, its flavour label as a superindex 
(and usually using Latin letters of the begining of the alphabet) and the "Clifford" label supressed; T is a 
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matrix that acts on Dirac indices (if it carries a Lorentz index we denote it with a Latin index of the middle 
of the alphabet). Powers of momenta (corresponding to derivatives in position space) are always indicated 
explicitly. For instance, the usual kinetic term would be 

We repeatedly use Fierz transformations in order to reduce bilinears to the diagonal in flavour indices form, 

S 12 =^^ 2 , P 12 = Vv i7s V£ 2 , V7 2 ^VV7^p 2 - (A.4) 

When we do so all flavour matrices disappear. 

Even more, for the sake of simplicity, we often drop the integral signs, the delta functions and the powers 
of (27r) d . However, when a term of the action is written they should always be understood to be present. 
Also, for the action in Appendix B it is convenient to define the sum and difference of the momenta of every 
bilinear. For instance, the above kinetic term will be written as 

\pi2^V( 2 (A.5) 

with p^ n = (p m ±p n y and with the integrals and momentum conservation delta function being assumed. 
The conventions for the Clifford algebra are the usual ones on the two-dimensional Euclidean space [39] , 

{ 7 \7 J } =2^', 7a = -» 7 V. (A.6) 

The completely antisymmetric tensor e y is e 12 = — e 21 = 1. 

We always use N to indicate the number of flavours and rj to denote the anomalous dimension of the 
fields. The conventions for the coupling constants are a Latin or Greek letter with a subindex. The letters 
are: 1) g for the four-fermions, non-derivative operators, 

5i - Pi 2 P 34 ) + toV&Vl ; (A.7) 

2) to, r, s and t for the various types of four-fermions, two-derivatives operator couplings; 3) a, c and e for 
the three types of six-fermions, one-derivative operator couplings; 4) Greek indices k, l, v, e and c for the 
six-fermions and three-derivatives operator couplings. 

When we expand the full action we introduce the scheme-dependent parameters 

Sometimes they enter solely in the combinations u = ^ , z = . 

A dot in any quantity means a derivative with respect to the RG flow parameter t. 
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Appendix B. The action 

In this appendix we present the complete action we use for the computation of the /3-functions. For the 
sake of clarity the action is divided into subactions according to the number of fermions and the number of 
derivatives, and in the case of six fermions and three derivatives also according to the fermionic structure. 
The integrals over the momenta and the (^-functions of global momentum conservation (with their respective 
powers of (27r) d ) are always omitted. The notation is presented in the Appendix A. 

s (2,1) = \v~JivL. 

S (4 ' 0) = 9i (S 12 S 34 - P 12 P 34 ) + toV&Vgt. 

S (4 ' 2) = {mip^ 2 + m 2 Pi2 ■ P34 + m 3Pi2 2 } x (S12S34 - P12-P34) + {rwX 2 + r 2Pi 2 ' P34 + r 3Pi 2 2 } x ^12^34 

+ WPuPli + S2Pi2°Pl2 k + S3Pl2 i P34 k + *4Pm Pu*} X + «PmW* (^12-P 3 4 - P12S34). 

= {a lV -i + a 2 p- 6 3 } x z(S 12 S 34 - PiMV^ + {<W + <W> * ^12^6 
+ e lP + V fe i(P 12 S 3 4 - S 12 P 3 4)Vi 6 . 

S (6 ' 3)a = {k!Pi 2 3 Pi 2 ■ p u + H 2 Pi2Pl2 ■ P34 + K 3Pi2 J P34 ' P 3 4 + ^Pu ?34 ' PfS6 + K 5Pl2 Pu ' P34 
+ ^Pl2P34 ■ P56 + K 7pt/Pl 2 ■ P34 + ^PHP1 2 ' P56 + «9?34 P34 ' P34 + P34, ' Pfti 

+ Kupt/Pu ■ pt 6 + K12P34 P56 ' P56 + ^PuPu ■ Pl2 + ^UPuPu ■ P34 + ^PuPu ' 2>56 
+ K16P34 J P34 ' P34 + K nP34 1 P34 ' P56 + ^ISPu P34 ' Ps~4 + ^WPu Ps~4 ■ p^ + K20P34 'P56 ' 2*56 
+ K 2lP34 5 P56 ' PfTe) X iVf2(S , 34S , 56 - fk^ki)- 

^(6,3)6 = { tl p- 2 Jp- 2 . p- + L2Pi2 3 Pl2 ■ P34 + '^2^34 ' P34 + ^^^34 ' P56 + L $P\2 P34 ' P34 
+ I-6P12P34 ■ P56 + wtiPvi ' P34 + wti ' Pl2 ' P56 + wtl P34 ' P34 + HoP^Pm ' P56 
+ tllP^Pii ' P56 + tl2P34 J P56 ' P56 + ^34 P\2 ' Pr2 + ^P34 P\2 ' P34 + HsPlwPw • p^ 

+ H6Pm P34 ' P34 + tirPw P34 ' Pm + *18Pm Pik ' Pik + *19Pm Pik ' Psa + Wj« ' Pse 

+ ^21^34^56 "Pie} X iVf 2 (S'34S , 56 ~ ^34^56)- 
5 (6 ' 3)C = {viPl2Pl2 k Pvt + v 2Pl 2 Pl 2 k Pl2 l + V3Pl2 i Pl2 k P34 l + t '^2 W W + v Wl2 P^Pls 
+ VeP^PlJ'PlJ + V7pthl2 k Pl2 + V S pi 6 j Pl2 k Pl2 + VvP34Pl2 k Pl2 + VwP^P^Pm 
+ WlW'PwW + ^Plhl^Pll + VY3ptip\2 k Pll + Vl4Pw 3 Pl2 k P34 l ) X ^^sV 

5 (6,3)<i = { eip +-7p+ . p + + eypviVt?. ■ P34 + £3Pi2 J P34 • P34 + £4Pr2 J Pr2 • Pr2 + £5Pr2 : 'pr2 • P34 

+ eeP^/P^ • P56 + £ 7Pr2 J P34 • P34 + £ 8Pr2 J P34 ' P56 + £ 9Pr2 J P56 ' P56 + ^lOP^P^i ' Pf2 
+ £llPr2 J Pl2 • P.34 + £l2Pr2 J Pl2 • P56 + £ 13Pi"2 J P34 ' Pr2 + ^UP^ P 34 ' P34 + ^UPm P34 ' P56 
+ £l6P56 J Pl2 ' Pf2 + ZnPh%P\2 ■ P34 + £ 18P56 J Pl2 ' Pbd) X (512^34 - P\2S 3 i)V^. 
5 (6 ' 3)E = {ftP^P^P^ + S2pi 2 3 Pl2 k P34 l + ^Pihl2 k P56 + W+jPuPto + fcPu'PuVa 
+ ^Pl2 3 P\2 k P34 + ^Pl2 3 Pl2 k P56 + ^P\2P34 k P\2 + ^Pl2 'PsV W + SwP\2 Psi Ptt 

+ ftlPse^W + Sl2P56 3 Pl2 k P34 l + CISPm PuW 1 X ^^(^12^34 - P^^)^. (B.l) 

The conventions to get rid of non-independent operators are as follows. In S^ 4,2 ) we consider terms 
containing multiplied by neither (S^S^ — P12-P34) nor V/2V34, because, for instance, ' P34^i2^34 = 
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— Pi2 2 ^i2^34- We do not take into account operators containing (S12S34 — Pi2P34)V^ 6 p^ 6 k , V k 2 V 34 V^ 6 p^ 6 l or 
(S12-P34 — Pi2»S'34)Vg' 6 p^ 6 ' c either. Finally, for terms like V( 2 V 34 we integrate by parts if they are multiplied by 
vti orp+ 4 \ e.g., V^V k 4 pf 2 j p+ 4 k = - V( 2 V k 4 p^p+ 4 k = V( 2 V k 4 pi 4 p+ k , and similarly for multiplied 

by Pi2 J , Pti OT PB6- 

We do not include S^ 2 ' 3 ) = iVj^CPi — P2) J '(pi ~~ P2) 2 as we have discussed in Sect. 3. Without it, we have 
a basis of 106 independent operators. 

Appendix C. (3 functions 

In this appendix we present the complete set of /3- functions. The scheme dependent parameters a, (3, 7 
and S have been defined in Eq. (5.3). 

■q = 4a[-mi + m 2 + m 3 + si + s 2 + s 3 + s 4 + t - 2N(r 2 + 2s 3 + s 4 )], 
gi = 2g 1 rj + 2a(oi + a 2 + ei - a 2 7V) + 2/?(-£i - S4- e 6 - e 9 - e w - e 12 - e 16 - £i 8 

+ 2Ki + K 2 + K 3 + K 5 + K 7 + Kg + K13 + K14 + K i6 + K i8 - 8K1./V), 

.92 = 2.92^7 + 2a(ai + ci + c 2 + e x - c\N - c 2 N) + (3[2(-ei + e 2 - £3 - £4 - £5 - £7 - £10 - £11 + £13+ 
ei 4 ) + 2vi + 2v 2 -v 3 + v 4 + v 5 + v 6 - v 7 - v$ + v 9 + 2(k 9 + k w - Ku - K12 + «16 - «17 + «18 

+ «19 + «20 + «2l) + 2(2ti + ii 3 + (44 + il6 + tl8 + «-2 + «-3 + &5 + «-7 + ^9) - 8iV(w 2 + 2ii + ti 8 )], 

rhi = (-2 + 2r\)m\ + ^[-2£i + 5£2 - 6e 3 + £ 6 + 2e 7 - 4e g + £i - 3ei 3 + £15 + e 16 + e u - 5ei 8 
+ ft + ft + 2^4 + ft - ft - fto - fti - ft2 - ft3 + 4ki - K 2 + 6k 3 - 6k 4 + 4k 7 - 5k 8 - Kg 
+ 2k 10 + Ku - K44 - «i7 + 2k 18 + 4k 20 + 4N(-2k 3 + 2k 4 - k 7 + k 8 )], 

OL 

m 2 = (-2 + 2?7)m 2 + ^[-£5 - £s + £11 + £17 - ft + ft + ft - 2 ft - 2ft + Ti2 
+ 3k 2 + 6k 6 — 3ki + 2«i 3 + 5k 15 + 3ki 9 — 2k 20 + 2k 2 i — 47V(2k 6 + K15)], 

OL 

m 3 = (-2 + 2r?)m 3 + ^ [~2£i - 4e 7 + £10 + £12 + 2ei 3 + e 16 + e 18 + ft + c 7 - 2q 8 + ? n + ci 3 

+ 4k 1 + 2k 2 + 6k 5 - k 7 — Kg + 2k 13 + 4k m + 2ki 8 + 2k 19 + Ak 21 — AN(2n§ + K14)], 
h = (-2 + 277)7-1 + - [2(-3£i - £ 2 + £ 3 + £ 4 + £5 - 3e 7 + e w + en + 3e 13 - e 14 ) 
+ 2v\ - 2v 2 — 3w 3 — v 4 + v 5 — 5v 6 - 3w 7 + 3t> 8 + v g + 4vu + 4w i3 

+ 2(-Kg + K10 - Ku + K12 + «16 + «17 + K 18 - «19 + «20 + «2l) + 2(4ii - (-2 + 6t 3 - 6i 4 

+ 4i 7 - 5t 8 - t 9 + 2tio + til - - ti7 + 2ti 8 + 4i 20 ) + 8iV(-2i 3 + 2i 4 - t 7 + t s - t 2 o)], 
r 2 = (-2 + 2r/)r 2 + — [-u 4 + w 6 - w 9 - 2w i2 - t>i 3 - 3wi 4 + 2(-/c 7 + k 8 + K14 + K15) 

+ 3i-2 + 6t6 — 3no + 2ii3 + 5ti5 + 3iig — 2t 2 o + 2i 2 i — 8N(i e + tis)], 
r 3 = (-2 + 2r7)r 3 + — [2(-3ei + 3e 2 - 3e 3 + £4 + £5 + £7 - 4e 9 + £10 + £11 - £13 - £14 - 4ei 8 ) + 2v\ - 2v 2 

-V3-V4- V 5 + V 6 + V 7 - V 8 - 5wg + 2(-Kg - K10 + Ku + K12 + 4ki 3 + Ki 6 - K17 + Ki 8 + Kig + K 20 

+ K21) + 2(4ii + 2t 2 + 6t 8 - t7 - tg + 2ti 3 + 4ti4 + 2n 8 + 2tiQ + 4t2i) - 8Ar(2i 5 + t 13 + <, 14 + t 2 i)], 
si = (-2 + 2t7)si + |[2(-£i + £3 + £4 - £7 + £11 + £13) - 2(?i + <r 2 + ft + ft) 
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+ 2i>i - 6v 2 - v 3 + v 4 + v b - 7v 6 - v 7 - v$ + vg + 8vu + 8w i3 

+ 2(-Ki - Kll + «16 - «18 - K 20 + «2l) + AN (v 6 - 2vu - V 13 + L W )], 

s 2 = (-2 + 2rj)s 2 + — [2(ei - £ 2 + £ 3 - £ 4 - £5 - £7 + 2£i 8 ) + 2(<r 5 + ft - ft - ft + 2?i 3 ) 

- 2vi + 6v 2 + v 3 + 7v 4 + 7v 5 - v 6 + v 7 + v 8 + 7v 9 

+ 2(2k 2 - Ki6 + K 17 + K 1S + «19 - «20 + «2l) - 4iV(w 4 + 2v 5 + Vg + L 2 + tig)], 
S3 = (-2 + 2r?)s 3 + |[2(-£ 6 - £ 8 + £12 - £15) + 2(-<r 3 - <T4 + ?7 - TlO + 2?n + 2ft 2 ) 
+ v 4 + Av 5 - v 6 + v 9 + 12w 10 + 4w n + 6v 12 + v 13 + 3v 14 

+ 4(-« 3 + k 4 + k 5 + k 6 ) - 4io - 2ii3 - <-15 + 4 19 + <-2 + 2i 6 + 2i 20 - 2i 2 i - 8iV(3wi + v 12 + L 6 )], 

OL 

s 4 = (-2 + 2i])s 4 + — [3(-l> 4 - v e + v 9 + 2v 12 + v 13 + 3v 14 ) + 2(-/t 7 + k$ + k 14 + k 15 ) 

- 42 - 2i 6 + tio + 2413 + 415 - 419 - 24 20 + 24 2i - 4iV(2wi 2 + 3^14 + 415)], 

OL 

t = (-2 + 277)4 + —(e 2 + 2e 3 - 3£ 5 - 3£ 6 - 2e 7 - 3e 8 - 4e 9 + 3en + 3£ i2 - £13 - 3£ i5 + 2£i 8 

+ 3(-?i - ft - ft - ft + ft + ft - ft - 2ft - fto + 2fti + 2ft 2 + 2ft 3 ) 

+ 2k 2 - 2k 3 + 2k 4 + 2k 5 + 2k 6 + k 7 — k 8 + k 10 — 3«n — 4ki 3 - «i 4 — k 15 + 3k 17 - k 19 — 2k 20 + 2k 2 i 
+ 4N(ea + e 8 - £12 + £15 + ft + <T4 - ft + fto - 2<m - 2ft 2 )), 

01 = (-2 + 3fj)oi + 4 ff i ff2 7, «2 = (-2 + 3r7)a 2 + 2^7, Ci = (-2 + 3r?)ci + 4.g| 7 , 
c 2 = (-2 + 3??)c 2 - 2 5 f 7 , e = (-2 + 377)e - 4^7 - 4.g l527 , 

Ki = (-4 + 3t?)ki + fy?/2 + 2.gi 7 m 3 , k 2 = (-4 + 37/)k 2 + 2 7ff i(2m 2 + s 2 ), 

k 3 = (-4 + 3t/)k 3 + <5.g? + 2 7ff i(2mi + m 3 + 4s 3 + i), k 4 = (-4 + 3t7)k 4 - Sg\ + 2-/ gi (2s 3 - t), 

k 5 = (-4 + 3?7)k5 + 2 7 5i(m 3 + 2s 3 ), k 6 = (-4 + 3r/)K e 

k 7 = (-4 + 3?7)K7 + 25g\ + 2j gi (2m 3 + 4r 2 + 4s 4 - t), 

k s = (-4 + 3r?)K 8 - 2Sgl + 2 7 gi(-2m 3 + 2r 2 + 2s 4 + t), 

kg = (-4 + 3^9 + 28gm + 2 7 (- 5 i(2si +t) + 2g 2 (m 2 + m 3 - t)), 
«io = (-4 + 3t])k w + 8S gi g 2 + 2"f(g 1 {2m 2 + 2r 3 - 2s\ + 3s 2 ) + 2.g 2 (m 2 + 2m 3 - t)), 
/til = (-4 + 3»?)kii + 4S gi g 2 + 2 7 (.g 1 (-2m 2 + 2r 3 - 2sj + t) + 2.g 2 (2m 2 + m 3 )), 

K12 = (-4 + 3?7)ki2 + 45gi.g 2 + 2 7 (#i(-2si + s 2 ) + 2g 2 (m 2 + 2m 3 - t)), k 13 = (-4 + 3r?)Ki 3 + 2 7 5ir 3 , 

k 14 = (-4 + 3??)ki 4 + 4 7 .gi(r 2 + s 4 ), k 15 = (-4 + 377)ki 5 , 

«16 = (-4 + 3t/)ki 6 + S 9l g 2 + 2 7 ( 5 i(2n + t) + g 2 (m 3 + 2t)), 

k 17 = (-4 + 3tj)k 17 + A5g 4 g 2 + 2 7 (.gi(4ri + s 2 - t) + Ag 2 {m 3 + t)), 

«18 = (-4 + 3?7)ki 8 + S gi g 2 + 2 7 (gi(r 3 + s 2 ) + g 2 m 3 ), 

«19 = (-4 + 3?7)ki9 + 4 7 g 2 m 2 , k 2 o = (-4 + 3r?)K 20 + 4<5gi.g 2 + 7 (fli(2ri + r 3 + 2s 2 ) + g 2 {2m 1 + m 3 + 2t)), 
«2i = (-4 + 3?])k 21 + 2-fg 2 m 3 , L x = (-4 + 377)41 - 5g 2 /2 - 2-fg 2 r 3 , L 2 = (-4 + 3r?)4 2 + 2 7 g 2 (-2r 2 + s 2 ), 
43 = (-4 + 377)43 - Sg\ + 2 752 (-2n - r 3 + 4s 3 ), 4 4 = (-4 + 37/)4 4 + Sg\ + 4 752 s 3 , 

23 



k 


= 


(-4 + 377 


k 


= 


(-4 + 377 


ho 


= 


(-4 + 377 


i-12 


= 


(-4 + 377 


hi 


= 


(-4 + 377 


hi 


= 


(-4 + 377 


^19 


= 


(-4 + 377 


421 


= 


(-4 + 377 


V3 


= 


(-4 + 377 


^6 


= 


(-4 + 377 


i>8 


= 


(-4 + 377 


Ull 


= 


(-4 + 377 


via 


= 


(-4 + 377 




= 


(-4 + 377 


£3 


= 


(-4 + 377 


£4 


= 


(-4 + 377 


£6 


= 


(-4 + 377 


£8 


= 


(-4 + 377 


£10 


= 


(-4 + 377 


£13 


= 


(-4 + 377 


£15 


= 


(-4 + 377 


£18 


= 


(-4 + 37/ 


<3 




(-4 + 377 


<7 




(-4 + 377 


iio 




(-4 + 37/ 


ii3 




(-4 + 377 



)i 5 + 2 192 {-r 3 + 2s 3 ), k = (-4 + 3fj)ta, h = (-4 + 3t])i 7 - 25g 2 2 + 4 7 .g 2 (2r 2 - r 3 + 2s 4 ), 
)t 8 + 2Sg 2 + 4jg 2 (r 2 + r 3 + s 4 ), i g = (-4 + 377)49 + 2Sg\ + 4 7 .g 2 (2r 2 + r 3 - Si), 
)tio + 85g% + 2 7ff2 (8r 3 - 2sj + 3s 2 ), ki = (-4 + 377)411 + A5g 2 + 4 7ff2 (3r 2 + 2r 3 - Si), 
)ti 2 + A5g\ + 2 7ff2 (2r 2 + 4r 3 - 2sj + s 2 ), k 3 = (-4 + 3i7)t 13 + 2 7 .g 2 r 3 , 
)ti4 + 4 7 .g 2 (r 2 + Si), L15 = (-4 + 3?y)ii5, he = (-4 + St])li 6 + Sg 2 + 2 7 g 2 (2ri + r 3 ), 
)ti 7 + (54.92 + 2 7 g 2 (4ri + 2r 3 + s 2 ), ks = (-4 + 3r/)t 18 + Sg 2 + 2 7 .g 2 (2r 3 + s 2 ), 
)tig + 4 7 .g 2 r 2 , i 20 = (-4 + 3r/)t 20 + 4.g 2 <5 + 4 7 .g 2 (2ri + r 3 + s 2 ), 
)t 2 i + 2 7 .g 2 r 3 , wi = (-4 + 3r/)?Ji - 4 7 g 2 s 3 , ?j 2 = (-4 + 3r])v 2 , 

)v 3 - 4 7 g 2 (si + 2s 3 ), 774 = (-4 + 377)^4 + 2 7 g 2 (s 2 - 2s 4 ), v 5 = (-4 + 3r/)7j 5 - 2jg 2 s 2 , 
)vq - 2~/g 2 (s 2 + 2s 4 ), v 7 = (-4 + 3t7)tj 7 + 4 7 .g 2 (si - 2s 4 ), 
)v s + 2~/g 2 (-s 2 + 2s 4 ), u 9 = (-4 + 3rj)v 9 + A^g 2 s 4 , v w = (-4 + 3r))v w , 
)vu - 2 7 # 2 (2si + s 2 ), V12 = (-4 + 3jj)ui2, t»i 3 = (-4 + 3t7)?ji 3 + 4 7 g 2 (si - s 4 ), 
)u u , ei = (-4 + 3t7) £i + S gi { gi + 7g 2 ) + 2j( gi (2 mi + m 3 + 2r x + 2r 3 ) + .g 2 (4m! + m 3 )), 
)e 2 + 28gi(-g 1 + 2g 2 ) + A-y{2g 1 ri + g 2 m 3 ), 

)£3 + - 2 52) + 27(51 (2mi + m 3 + 2n - r 3 ) + g 2 (-2m 1 + m 3 )), 

)e 4 - <5.gi.g 2 + 2j(g 1 (m 3 - r 3 ) + 5 2 m 3 ), e 5 = (-4 + 37/)e 5 + 4jg 2 m 2 , 
)e 6 + 4 7 gi(m 2 - r 2 ), e 7 = (-4 + 3r/)e 7 + 2Sg!g 2 + 2 7 (gi(m 3 + 2r 3 ) + g 2 m 3 ), 
)e 8 + 4 7ff i(m 2 + 2r 2 ), £9 = (-4 + 3r/)e 9 + 5g\ + 2 7 .gi(2m 3 + r 3 ), 
)eio - 2(5.gig 2 - 4 7 5 2 m 3 , £11 = (-4 + 3r/)en - 4 7 g 2 m 2 , e i2 = (-4 + 3?7)ei2, 
)ei 3 - 45gig 2 - ^(gm + g 2 m 3 ), i u = (-4 + i-q)e u - 4-fg 2 m 2 , 
)ei 5 - 4 7 gir 2 , £i 6 = (-4 + 3?7)ei 6 , ei 7 = (-4 + 377)ei 7 + 475im 2 , 
)ei 8 + 2Sg 2 1 + 4 7 .g 1 m 3 , ft - (-4 + 3r?)<ri + 2 7 ( 5 i(4si + s 2 ) + 3g 2 t), <t 2 = (-4 + 3r?)<r 2 , 
)?3j <T4 = (-4 + 3t7)<T4 - 4 7 .gis 4 , <r 5 = (-4 + 377)^5 + 2 7 gis 2 , <r 6 = (-4 + 3j?)? 6 , 
)<T 7 + 4 7 5is 4 , = (-4 + 377)^ 8 + 2 7 (.gis 2 -,g 2 i), ft = (-4 + 3r/)^ 9 - 2 7 g 2 t, 
K10 + 8 7 gis 4 , ?n = (-4 + 3?j)?ii +2 7 gi(2s 3 +t), ft 2 = (-4 + 3t7)<ti 2 - 4 7 gis 3 , 
) ft3 + 2 7ffl (- S2 + i). (C.l) 



The only manifestation of the cut-off function comes through the constants a, f3, 7 and S. However, as 
it happens in the bosonic case [19] we can reduce the number of independent parameters from four to two. 
In fact, by performing the following rescalings 



,(4,0) 



1 

a-f~ 



. fl (4,0) j 5 (4,2) 



07 



2 5 (4 < 2) , 5 (M) 



-^5 (M) , 5 (6 < 3) 
a z 7 



-2~2 9 



(6,3) 



(C.2) 



where by 

g (m,n) we 

denote the coupling constants corresponding to the operators with m fcrmions and n 
derivatives, it can be shown that the /3-functions depend on the scheme only through the combinations 



z = 4j and w = Moreover, z enters the equations only as a global factor of the anomalous dimension. 
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Figure Captions 

Fig. 1. Ni] (solid line) and Ai (dashed line) as functions of N. This solution matches with the Type I 
solution of the large N limit. 

Fig. 2. Ai as a function of w, (z = 0.5), for N — 3 and N — 10. The minimum clearly decreases with N. 

Fig. 3. rj (solid line) and Ai (dashed line) as a function of TV for z = 0.5 and w = —2. In both curves the 
upper branch corresponds to the solution that matches with the Type II large N solution. 

Fig. 4. 77 (solid line) and Ai (dashed line) as a function of N, (z = 0.5, w — —2) for a different fixed point 
solution. In this case both exponents are of order 1 as N — > oo. 
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